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Department  of  Physics,  University  of  Houston,  Houston,  Texas  77204-5506,  U.S.A. 

ABSTRACT 

We  propose  a generalized  ellipsometric  technique  using  a rotating  sample.  The  ellipsometer 
consists  of  a polarizer,  a rotatable  sample  holder,  an  analyzer,  and  a detector.  Fourier 
coefficients  are  measured  and  used  to  extract  the  system’s  dielectric  tensors  and  film  thicknesses. 
The  main  advantage  of  the  technique  is  that  all  parts  of  the  ellipsometer  are  fixed  except  the 
sample,  whose  azimuth  angle  can  be  modulated.  We  show  calculated  responses  to  isotropic  and 
anisotropic  materials  as  well  as  superlattices.  Potential  applications  for  characterizations  of 
anisotropic  nanostructures  are  discussed. 


I.  INTRODUCTION 

In  a conventional  rotating  element  ellipsometer,  only  one  variable,  i.e.,  the  azimuth  angle 
of  the  rotating  element,  is  changed  in  the  characterization  of  an  isotropic  material.  In  current 
generalized  ellipsometers  for  characterizing  anisotropic  systems,  however,  either  additional 
variables  [1],  e.g.,  the  angle  of  incidence  and  the  azimuth  angles  of  the  polarizer  (or  analyzer) 
and  of  the  sample,  are  changed,  or  components  not  commonly  found  in  a conventional  rotating 
element  ellipsometer  are  needed  [2,3],  Even  in  the  simplest  case  of  uniaxial  materials,  two 
variables,  i.e.,  the  azimuth  angles  of  the  polarizer  and  analyzer,  are  required  if  sets  of  V and  A 
are  measured  to  determine  the  dielectric  tensors  of  the  materials  [4].  Optimization  of 
ellipsometric  setups  with  fewer  variables  involved  in  measurements  is  highly  desirable  in 
situations  where  measurements  are  remotely  controlled,  since  the  reliability  of  the  controls  and 
the  accuracy  of  the  measured  data  are  improved.  One  such  situation  is  the  characterization  of 
samples  grown  in  space.  In  this  paper  we  present  a theoretical  development  and  show  that 
changing  only  the  sample’s  azimuth  angle  is  sufficient  to  determine  dielectric  tensors  and  film 
thicknesses  of  arbitrarily  anisotropic  systems.  Instead  of  *P  and  A,  the  intensity  dependence  on 
the  sample’s  azimuth  angle  is  measured.  Two  approaches  for  analysis  are  proposed  to  determine 
the  dielectric  tensors  and  film  thicknesses  from  the  measured  intensities.  Popular  conventional 
ellipsometer  setups,  e.g.,  polarizer- compensator- sample -analyzer  (PCSA)  or  polarizer-sample- 
analyzer  (PSA),  can  use  this  approach  to  characterize  anisotropic  systems  by  keeping  all 
components  fixed  except  the  sample. 


II.  THEORY 

In  this  section,  we  study  the  dependence  of  intensity  on  as,  die  sample’s  azimuth  angle, 
to  permit  extraction  of  the  dielectric  tensors  and  film  thicknesses.  For  an  n-layer  anisotropic 
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system,  following  the  same  procedure  and  coordinates  setup  as  in  Refs.  [4,5],  the  Jones  matrix 
( r ) for  reflection  ellipsomctry  is 


={T„T»-T32T2l)/{Tnr21-TnT2]), 

(1) 

= {TyiT»-T3iTX2)l{TnTu-TuT2,), 

(2) 

= {T4IT22-T42T2>)/(TUT22-TI2T2I), 

(3) 

= ( TA2Tn  ~ TJn )l{TuTn  ~ TnTn)  - 

(4) 

where  (T)  is  the  transfer  matrix  of  the  system.  The  formula  for  calculating  (7)  can  be  found 
elsewhere  [6].  These  expressions  establish  the  relationships  between  the  Jones  matrix  and  the 
system’s  dielectric  tensors  and  film  thicknesses.  For  layer  /',  if  we  use  the  diagonal  tensor  (£/),, 

/=  1,  2,  3,  to  represent  the  dielectric  tensor  with  respect  to  its  principal  axes.  0,=(0,„  0,.,-,  <p,„  0,,-, 

0W , (pii)  for  the  Euler  angles  of  the  frame  of  principal  axes  with  respect  to  the  laboratory  frame  for 
the  real  and  imaginary  parts  of  the  dielectric  tensor,  and  d,  for  the  thickness,  ( T)  is  completely 
determined  by  (£/);,  0/,  and  d,-,  (/'=  I , . . . «). 

If  the  PSA  setup  is  used,  the  intensity  after  the  analyzer  is 

=/J(r„,cosaP+t;;,sinaP)cosaA  +(r;„  cos «,,+/;  sinar)sinaA|  , (5) 

where  Id  is  the  intensity  reaching  the  detector,  Ir  is  a constant,  and  aP  and  aA  are  the  azimuth 
angles  of  the  polarizer  and  analyzer,  respectively. 

If  ar  and  aA  are  fixed,  changing  as.  i.e.,  rotating  the  sample,  leads  to  the  change  of 
0,  = (0,,,0[,).  Recall  that  the  transfer  matrix  (7)  depends  on  0, . Therefore,  the  Jones  matrix 
changes  with  as  and  so  docs  the  intensity.  However,  the  dielectric  tensors  on  the  principal  axis 
frame  and  the  Euler  angles  0,  and  0,  remain  the  same.  Substitution  of  Eqs.  (1)  - (4)  into  Eq.  (5) 
gives  Id  a function  h of  the  dielectric  tensors  and  film  thicknesses.  That  is, 

7rf(as)=  I,M(ei)h  0,-,0+as,  ©y,  0//,  0,/o+  as,  ©/,  dj),  (6) 

where  the  underlined  parameters  do  not  change  when  the  sample  rotates,  and  0,y<;  and  0„o  arc  the 
initial  angles  for  0,y  and  0,,-. 

From  Eq.  (6),  if  the  total  number  of  (£/),-.  0,  /,  0,  ,o,  <p,  i.  On-  0//o.  <pa . and  c/,  is  k,  and  7;,  is 
taken  as  unknown,  they  can,  in  principle,  be  determined  by  k+ 1 independent  equations  between 
Id  and  Ip,  (£/),-,  0,  ;,  0,-/0,  0,-/,  0,„  0,70,  0,;,  and  d,.  In  the  rest  of  this  section,  W'e  discuss  two 
approaches  of  changing  as  to  establish  m (m  > A’  + l ) equations  between  ld  and  /;„  (£/),.  0,„  0„o, 

0,/.  07,  0/70,  0/7,  and  dp 
Approach  one'. 

The  intuitive  way  to  establish  m (m>k  + 1)  equations  between  /,/  and  Ip,  ( £/),- , 0,„  0,,o, 
0,/,  0,/,  0,,o,  0,„  and  d,  is  t)  set  as  to  m different  settings.  If  the  intensities  7,/(asi),  1, /(as?),  ..., 
7,/(as,„)  have  been  measured,  m equations  are  set  up.  They  arc 

7;, //((£/)/,  gi , 0,/o+asi , 0,-/,  0,/,  0,-n+as i , 0//,  d/f  7,/(aSi ),  (7) 
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ipHtefih  Sn , A*+a&,  8th  ifeo+a©,  &h  di)=  4riasz),  (8) 


Iph((£})h  S'h  &}rii  8fi>  *2tf>  ^i)-  (9) 

To  solve  these  equations  for  (£/),-,  0,„  <pri0,  (prl,  On , </>m , <pu,  and  d,-,  a test  function  is 
introduced: 


^ = X[/p*(fe,)p©„.f,o+«S,TP„^n>^o+^/P,,.<)-/I/K)]VCTj2’  (10) 

;=i 

where  <r;  is  the  standard  deviation  for  Id{a^).  Ideally,  if  there  are  no  errors  in  the 
measurements  of  I,i,  Eq.  (10)  reduces  to  Eqs.  (7)  - (9)  if  Err  is  minimized  to  zero.  In  practice, 
errors  exist  in  Id,  but  Err  can  be  reduced  to  a minimum  value.  The  Ip,  ( £/),- , 0,„  0,.,o,  <p >„■,  0<„  ipuo,  <Pu , 
and  di  that  minimize  Err  can  be  taken  as  the  true  values.  We  outline  the  algorithm  for  the 
minimization  and  discuss  the  error  distribution  in  determining  (£/),-,  0„  0„o,  (p,,,  da,  <pm , <Pu,  and 
di. 

For  clarity  of  notation,  we  use  xi,  I = 1,  2,  ...,  k,  k+\  to  represent  the  k variables  of  (£/),-, 
dri,  <j>,io,  <p,i,  6a,  Quo,  (pi i,  d,,  and  Ip,  and  define  a vector  XT=(xt,  X2, ....  Xjt+i),  where  the  superscript  T 
denotes  the  transpose  of  the  column  vector  X When  Err  reaches  its  minimum, 

(dErr/dxl)i  = 0,1=  1,2, ...,  fc+1,  (11) 


where  the  subscript  I means  that  xi,  i.e.,  (£/)/,  0„,  (f>,i0,  <pn,  6a,  (f>M,  <j%,  and  dt,  assume  the  true 
values  of  the  dielectric  tensors  and  film  thicknesses.  Using  matrix  notation,  Eq.  (11)  can  be 
simply  denoted  as 

VErr  = 2/4rH  = 0 , (12) 


where 

h =(AA> -A)r.  (13) 

hj  = IJt({,£,)l,0ri  ,0riO  + aSJ  ,<p  0„ , <pli0  + ,<pu,d,)-Id(aSJ),j=\,2,...,m,  (14) 

AJt  =(dhJ/dxl)J<Jj,j=l,2,...,m;  l=\,2,...,k+\.  (15) 

Near  the  minimum  of  Err,  H can  be  linearized  as  H=  H,  +AAX , and  substituting  it  into  Eq.  (12) 
obtains 

AX  = -(ATAy' ArH,=-MATH,.  (16) 


Eq.  (16)  gives  an  expression  for  changing  X recursively  to  minimize  Err.  The  standard  deviation 
of  xi,  1=1,  2,  ...,  k,  k+\,  determined  with  the  above  algorithm,  is  [8] 
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(17) 


< = M„. 

Approach  two: 

Instead  of  setting  a sample  to  different  discrete  angles,  the  sample  can  be  rotated  at  a 
constant  frequency  CO.  That  is,  as  = COt . In  this  case,  Eq.  (6)  can  be  expressed  as 


cos kot+h,  sin  tar). 


(18) 


where  a / and  hi  are  the  Fourier  coefficients. 

If  m=2j+\,  m>k  + 1 , Fourier  coefficients  oq,  a i,  b\,  ...,  «/,  /»,,  are  measured,  equations  of  (f/)„ 
Bri,  tfiriOj  (Pit*  On,  , (jO,,-,  and  can  be  established.  They  are 


jr  +«■  0„,<p„ .„  +a,  </>„,</,  )rto  = n0, 

Ui> 

yr  jlh  ((£/)/  ,0ri,0r,o  +a,?>„,e,i.^o+a.«P11.</i)cos  ./«</«  = a, . 

tJp 

( (£( ),> ®r, . 0a .,  + a > <P,,> + a, <)>,, , r/,) S i n yaz/a  = b, , 


where  f < 2;r  is  the  smallest  period  of  function  h. 
To  solve  Eqs.  (19)  - (21),  a test  function 


Err-t 


2 f h cos  lada 

Jo 


(1 + SJII, 


fa +2 


2 J /z  si  n lot  dot 


(19) 

(20) 
(21) 


(22) 


is  introduced,  where  cr/  are  the  standard  deviations  for  the  Fourier  coefficients. 

The  algorithm  outlined  in  Approach  one  can  be  employed  to  solve  Eqs.  (19)  - (21),  if  the 
intensities  are  replaced  by  the  Fourier  coefficients  and  the  function  h replaced  by  the  Fourier 
integrals.  Eq.  (17)  can  also  be  used  to  estimate  the  error  distributions  of  (£/)„  0,„  <pr, o,  t p„,  0,„  <p„o, 
(pa,  and  cli. 

If  j is  large,  direct  measurements  of  the  Fourier  coefficients  may  be  inconvenient.  In  this 
situation,  the  intensity  dependence  on  the  sample’s  azimuth  angle  can  be  measured.  With  this 
relation,  the  Fourier  coefficients  can  be  obtained  numerically.  It  should  be  noted  that  if  t< In,  all 
as  selected  in  Approach  one  should  not  be  different  by  t.  If  t= 0,  i.e.,  /,/  is  constant,  only  «o  is 
meaningful.  In  this  case,  this  method  fails. 


III.  MODELS  CALCULATION 

In  this  section,  we  prescribe  the  dielectric  tensors  and  film  thicknesses  of  some  common 
structures  and  conduct  forward  calculation  for  the  intensity  dependence  on  sample’s  azimuth 
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angle.  The  purpose  of  the  study  is  to  show  that  different  structures  possess  different  relationships 
between  intensity  and  sample’s  azimuth  angle.  Therefore,  from  a measured  intensity  dependence 
on  sample’s  azimuth  angle,  a backward  calculation  can  extract  the  dielectric  tensors  and  film 
thicknesses  by  use  of  the  two  approaches  described  in  the  previous  section. 

In  the  study,  we  assume  that  the  angle  of  incidence  is  70° , and  the  wavelength  is  600nm.  All 
materials  are  transparent  for  simplicity.  A PSA  setup  is  used  and  a?  = CtA  = 45° . An  isotropic 
material,  two  bulk  anisotropic  materials  (uniaxial  and  biaxial)  are  studied.  The  isotropic  and 
uniaxial  materials  are  also  used  as  the  substrates  of  a superlattice  structure,  which  is  20  periods 
of  a two-layer  structure  on  a substrate.  They  can  be  denoted  as  ABAB. . .ABC,  where  A and  B are 
biaxial  films,  and  C is  substrate.  The  properties  of  these  materials  are  listed  in  Table  I.  Relations 
between  intensities  and  sample’s  azimuth  are  shown  in  Figure  1 . Parts  of  the  Fourier  coefficients 
are  listed  in  Tables  II. 

Figure  1 shows  that  the  intensity  relationships  are  different  for  different  structures.  Table  II 
shows  that  except  for  isotropic  materials,  all  other  structures  contain  more  non- zero  Fourier 
coefficients  than  unknowns.  Therefore,  Eq.  (22)  can  be  used.  For  isotropic  materials.  Figure  1 
shows  that  the  intensity  is  a constant.  As  stated  in  the  above  section,  this  method  fails  in  this 
case.  One  interesting  observation  is  that  in  Figure  1,  two  superlattices  show  difference  mainly 
near  extremum  positions.  However,  Table  II  shows  that  all  Fourier  coefficients  have  distinct 
difference.  This  fact  implies  that  if  Approach  one  is  used,  as  has  to  be  set  near  extremum 
positions  so  that  the  difference  can  be  detected.  However,  Approach  two  does  not  have  this 
limitation. 


Table  I.  Refractive  indices  of  different  materials  along  their  principal  axes,  Euler  angles  of  the 
principal  axes  relative  to  the  laboratory  frame,  thicknesses  of  the  materials  and  number  of 
parameters  that  can  be  determined. 
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Table  II.  Fourier  coefficients  a;  and  bi  for  the  intensity  7 , (ft»)  = I cos  tot+6,  sin  tot). 


Lattice  1 (resp.  2)  is  a 40- layer  periodic  superlatticc  on  an  isotropic  (resp.  uniaxial)  substrate. 
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Figure  1.  Intensity  dependence  on  sample’s  azimuth  angle  for  ( I ) an  isotropic  material,  (2)  a 
uniaxial  material,  (3)  a biaxial  material,  (4)  a 40- layer  periodic  superlattice  on  the  isotropic 
substrate,  and  (5)  a 40- layer  periodic  superlattice  on  the  uniaxial  substrate. 

IV.  SUMARY 

In  this  paper,  we  present  a theoretical  development  to  optimize  a generalized  ellipsometer 
so  that  only  the  sample’s  azimuth  angle  needs  to  be  changed  in  the  determination  of  the  dielectric 
tensors  and  film  thicknesses  of  a arbitrarily  anisotropic  systems.  Five  models  arc  calculated  as 
examples. 
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